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We consider a class of multi-component hybrid inflation models whose evolution may be analyt- 
ically solved under the slow-roll approximation. We call it multi-brid inflation (or n-brid inflation 
where n stands for the number of inflaton fields). As an explicit example, we consider a two-brid 
inflation model, in which the inflaton potentials are of exponential type and a waterfall field that 
terminates inflation has the standard quartic potential with two minima. Using the SN formalism, 
we derive an expression for the curvature perturbation valid to full nonlinear order. Then we give an 
explicit expression for the curvature perturbation to second order in the inflaton perturbation. We 
find that the final form of the curvature perturbation depends crucially on how the inflation ends. 
Using this expression, we present closed analytical expressions for the spectrum of the curvature 
perturbation Vs(k), the spectral index ns, the tensor to scalar ratio r, and the non-Gaussian pa- 
rameter /]v£ a \ in terms of the model parameters. We find that a wide range of the parameter space 
(ns, r, /jv£ al ) can be covered by varying the model parameters within a physically reasonable range. 
In particular, for plausible values of the model parameters, we may have a large non-Gaussianity 
/jvzf 1 ~ 10-100. This is in sharp contrast to the case of single-field hybrid inflation in which these 
parameters are tightly constrained. 

PACS numbers: 98.80.-k, 98.80.Cq 

6 

t*. 

The standard, single-field slow-roll inflation predicts the curvature perturbations which are almost Gaussian to high 
accuracy [1, 2], Thus the detection or non-detection of non-Gaussianity will have extremely important implications 
to theories of the early universe, and a variety of multi-component models that give detectable non-Gaussianities 
have been proposed [3-10]. In many of these cases studied so far, however, it is not quantitatively clear how the 
non-Gaussianity arises and what determines its level. 

In this paper, we study a class of multi-component hybrid-type inflation models whose dynamics can be analytically 
solved. Similar exactly soluble slow-roll models were investigated in Ref. [11]. Here we extend the analysis given in 
Ref. [11] by including a coupling to a water-fall field which terminates the inflationary stage. Then applying the SN 
1^ | formalism [12-15], we compute the curvature perturbation analytically. Models of hybrid inflation similar to the one 
presented in this paper were proposed and investigated by Bernardeau and Uzan [5] and by Alabidi and Lyth [9] . 

In the <WV-formalism, the final amplitude of the curvature perturbation on comoving hypersurfaces 1Z C (or equiv- 
alently on uniform total density hypersurfaces £) is given by SN, where SN is the perturbation of the number of 
e-folds between the initial flat time-slice at horizon crossing during inflation and a comoving time-slice during the 
final radiation-dominated stage by which time all the inflationary trajectories are assumed to have converged to a 
unique one [13, 15]. 
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I. INTRODUCTION 



For the models we consider in this paper, assuming the scalar field perturbations are Gaussian up to the time of 
horizon crossing, the non-Gaussianity is local in the sense that the curvature perturbation at each spatial point is 
expressed in terms of a nonlinear function of Gaussian fields at the same spatial point. In this case, the level of 
non-Gaussianity is conveniently characterized by the quantity /]^£ al [16]. 

This paper is organized as follows. In Section II, we first consider conditions for a slow-roll model to be exactly 
soluble and derive a general formula for the curvature perturbation. In Section III, as a specific example, we present 
a model of multi-field hybrid inflation, which we call "multi-brid inflation". Then focusing on the case of two-brid 
inflation, we give various formulas explicitly. In Section IV, we compute analytically the spectrum of the curvature 
perturbation Vs{k), its spectral index ns, the tensor-to-scalar ratio r, and the non-Gaussian parameter /j^£ al . They 
are all expressed explicitly in terms of the model parameters. We conclude the paper in Section V. Some computational 
details are described in Appendix A. We use the Planck units where = 8nG = 1. 
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II. EXACT SOLUBLE CLASS 



We consider an Einstein-scalar Lagrangian, 

l = \r- \g^h ab {^)d^ a dA b - V{4>) , (2.1) 

where R is the Ricci scalar, the Latin indices a, b run from 1 to M, h ab is the field space metric, and V is the potential. 
We assume the dynamics of inflation is described by this Lagrangian. Later we introduce a water-fall field x that 
terminates inflation. For the moment, however, we concentrate our discussion on the inflationary stage. 
The Friedmann equation and the field equations are, respectively, 

3H 2 = l -h ab ^ a ft + Vtf) , 

4> a + 3H<j) a + h ab d b V = , (2.2) 

where H = a/a with a being the cosmic scale factor, and a dot denotes a derivative with respect to the cosmic proper 
time t, ' = d/dt. The slow-roll equations of motion are obtained by neglecting the kinetic term in the Friedmann 
equation and the second time derivative in the field equation, 

3H 2 = V{cj>) , 

3H<j) a + h ab d b V = . (2.3) 

For later convenience, we change the time variable from t to the number of e-folds from the end of inflation backward 
in time, 

dN = -Hdt . (2.4) 

Then the slow-roll equations of motion give 

d<p a h ab d b V h ab d b V 



dN 3H 2 V 



(2.5) 



Now we consider the case when the slow-roll equations of motion (2.5) can be exactly soluble. They can be exactly 
solved if the right-hand side of Eq. (2.5) for each index a takes the form, 

h ab d b v _ f(f) i0 

where F represents an arbitrary function of (<^ 1 , <j> 2 , ■ ■ ■ , <p M ), and f a is a function of only <fr a for each a (a = 
1, 2, • • • , M). In this case, we have 

1 ,7J,a 1 

(2.7) 



f a (<P a ) dN F{<t>) ' 
Introducing a new set of field coordinates q a by 

>■«•-/£. m 

the equations of motion become 

*a±!L = L. (2.9) 

dN F y ' 

Here we introduce the radial and angular coordinates of the field space, 

q a = qn a ; ^(n») 2 = 1 . (2.10) 

a 

Then it is straightforward to see that n a is conserved, 
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FIG. 1: A schematic diagram of classical trajectories in the field space with the coordinates q a . The angular coordinates 
n " = <f ll are conserved, and hence all the trajectories are radial in these coordinates. The curve indicated by q — qj denotes 
the surface at which the inflation ends, which may depends on n a . The e-folding number N is counted backward in time from 
Q = If- 

and Eq. (2.9) reduces to a single equation for q, 

ding 
dN ~ 

where F is now regarded as a function of q and (n , n 2 , 



1 



F(q,n a ) ' 



(2.12) 



In fact, what we have shown is not restricted to slow-roll inflation. Whenever a system approaches an attractor-like 
asymptotic stage, the value of the field at one instant of time completely determines the motion of the field at its 
subsequent time. In such a case, we have M first-order differential equations instead of M second-order differential 
equations. Then the exact solubility means there are M — 1 constants of integration. This is indeed the case we have 
considered in the above (there are M — 1 degrees of freedom in n a ). A visualization of this situation is depicted in 
Fig. 1. 

At this point, let us consider a couple of somewhat more specific examples. The condition for exact solvability (2.6) 
is satisfied if one can choose a set of field space coordinates such that the field space metric takes the form, 

,q fc _ diag. (h'm,hH^),---,h-(r)) 

H(<t>) ' ( j 



and the potential is ether of product type, 



in which case we have 



V = l[V a (^), (2.14) 



d\u V a 

f 1 = h a , F = H. (2.15) 



or of separable type, 



a 



(2.16) 
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in which case we have 



dV a 

f a = h a —, F = HV. (2.17) 



The separable case was first investigated in detail by Starobinsky [17]. 

Now, going back to the general case, Eq. (2.11) implies the trajectories are always radial in the field space spanned 
by q a , and Eq. (2.12) can be easily solved for N to give 

N(q,n a )= f Q F{q',n a )d\nq' , (2.18) 

Jq f 

where qf is the value of q at the end of inflation. Note that qf generally depends on n a . Thus, applying the nonlinear 
SN formula, we immediately obtain a fully nonlinear expression for the conserved comoving curvature perturbation, 

SN = N(q + 5q, n a + 5n a ) - N{q, n a ) . (2.19) 

It is instructive to consider the linear limit of the above SN. Denoting it by 8lN, we find 

S L N = F {q ,n^+ f 9 ^^.LpL Sn a i 
q J q on a q' qj dn a 

where the values of (q, n a ) are those at horizon crossing during inflation and (5q, 8n a ) are the fluctuations evaluated 
on the flat hypersurface at that epoch. It may be noted that 5q a are given by the fluctuations of the original field 
variables 5(f) a as 

S6 a 

S\nq a = -^—. (2.21) 

The formula (2.20) consists of three terms. The first term may be interpreted as the one due to adiabatic perturbations, 
the second to entropy perturbations during inflation [13, 18, 19], and the third to entropy perturbations at the end of 
inflation [6-8]. It may be worthwhile to mention that although these distinctions are meaningful and useful at linear 
order (or perhaps perturbatively), they are not so at full nonlinear order. 

Figure 2 is a schematic graph to explain the meaning of these three terms. The thick lines represent three different 
kinds of trajectories in the field space, and the thin lines with arrows on both ends represent the field fluctuations. 
The wavy dashed line represents the surface at which inflation ends. After inflation, each kind of orbit converges to 
a unique trajectory. 



The one on the left, where the fluctuations are parallel to the orbits, corresponds to the first term in Eq. (2.20). 
In this case, the fluctuations in the initial condition directly gives SN. This is identical to the conventional adiabatic 
curvature perturbation in the single field case. The one in the center corresponds to the second term. The fluctuations 
are orthogonal to the orbits, so they are entropy (isocurvature) perturbations during inflation. However, by the time 
the two adjacent orbits converge to a unique orbit, one finds the number of e-folds depends on which route the universe 
has taken. This gives rise to SN before the end of inflation. The one on the right corresponds to the third term. The 
entropy perturbations are not converted to SN until the end of inflation. However, the surface that determines the 
end of inflation may not be orthogonal to the orbits. This gives rise to SN in the end. 



III. MULTI-BRID INFLATION 



We consider a multi-component inflaton field with the potential, 

M 



V(0) = y o exp 



A=l 



(3.1) 



Without loss of generality, we assume rriA^A > for all A = 1, 2, • • • , M. The slow-roll equations of motion are 

d(j) A 1 dV 



dN V dd> A 



m A ■ 



(3.2) 
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FIG. 2: Schematic graphs describing the three terms in Eq. (2.20). The thick lines represent three different kinds of orbits in 
the field space, and the thin lines with arrows on both ends represent the field fluctuations. The wavy dashed line represents 
the end of inflation. The one on the left corresponds to the first term, the one in the center to the second term, and the one 
on the right to the third term. See text for more details. 



Note that the effective mass-squared M\ for each <f) A is given by 



d 2 v 



= m A V = 3m 2 A H 2 . 



Thus the slow-roll condition is satisfied for m A <C 1. Introducing a new variable q A by 

<L4 = cxp[<j) A /m A ] = qn A ; ^(n A ) 2 = I , 

A 

the slow-roll equations become 



d In q ^ dnA _ Q 



dN 



dN 



These are immediately integrated to give 



N = In q — In qj = — In 



\ ' e 2<p A /m A 
. A 



2 



/m A 



Now, similar to the conventional single-field hybrid inflation [20], we assume that inflation ends at 

^9a4>a = ° 2 i 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



and the universe is thermalized instantaneously. This may be realized by introducing a water-fall field \ that termi- 
nates inflation. Specifically, we promote Vq in Eq. (3.1) to a function of x (and of <p A ) as 



Vo = IJ2 9a Kx 2 + j(x 2 -^ 



(3.8) 



Then we see that the field \ IS stable at the origin for 'Y^ A g A $ A > & , but becomes unstable for J^aSa^a < a 
which brings inflation to an abrupt end. 
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Now we restrict our discussion to the case of two inflaton fields, that is, two-brid inflation. Then it is convenient to 
parametrize (A = 1,2) as 



n\ = cos 9 , n 2 = sin# . 
Similarly, we parametrize the scalar field at the end of inflation as 

<Pi,/ = — cos 7 , 02,/ = — sin 7 . 
9i 92 



(3.9) 



(3.10) 



Note that since rriA&A > (A = 1,2), both 77i 1 g 1 cos7 and TO 2 g2sin7 are positive. Then, since 9 is a constant of 
motion, we have 



In 



a cos 7 cr sin 7 



mi rri2 gimi 



92m 2 



(3.11) 



This equation may be solved for 7 in terms of 0i and 02- Hence we may regard 7 as a function of (0i, 02), 7 = 7(0i, 02)- 
Therefore, the number of e-folds (3.6) may now be regarded as a function of (0i, 02), 



JV = JV(0i,02) = 2 ln 



o 20i/m 



1 +e- 



202/m2 



g 2(7 cos 7/(91 mi) _|_ e 2ersin7/(g 2 m 2 ) 



It is then straightforward to obtain <5A to full nonlinear order. It is given simply by 

SN = iV(0i + (50i, 02 + <50 2 ) - AT(0i, 2 ) . 



(3.12) 



(3.13) 



Here it may be useful to clarify the origin of SN in the present model. Going back to the general expression for the 
linear curvature perturbation (2.20), and noting that we have F = 1 in the equation of motion (2.12) in the present 
model, we see that there is no contribution from the second term which is due to entropy perturbations during 
inflation, Thus the linear curvature perturbation is due to the initial adiabatic perturbation at the time of horizon 
crossing and to the final entropy perturbation at the end of inflation. Explicit computations of these contributions to 
SN are given in Appendix ??. 

Before closing this section, we note that there is a correction to be added to the above formula in the rigorous 
sense. It comes from the fact that the surface where the inflation ends in the field space, Eq. (3.7), is not equal to 
the surface of constant energy density. This means that the hot Fricdmann stage of the universe starts at slightly 
different temperatures for different values of (0i./, 02,/) at which inflation was terminated. This can be taken care of 
by introducing a correction to the number of e-folds given by Eq. (3.12). Namely, denoting the potential energy at 
the end of inflation by Vf, we add a correction term N c to Eq. (3.12) as 



N= - In 
2 



o 20i/mi _j_ g20 2 /TO2 



e 2<r cos 7/(91 mi) _|_ e 2<jsin7/(g 2 m 2 ) 



+ AT C (7) ■ 



where 



N c = - In 
4 



'Vf' 




To. 





mi fn-i . 

— cos 7 H sin 7 

9i 92 



(3.14) 



(3.15) 



In the above we have assumed that the universe has become radiation-dominated right after inflation. It is then easy 
to see that this correction term is negligible for sufficiently small mi and mi. To be more specific, if we compare the 
7-dependent terms in the original N to those in N c , the dependence is apparently strong in the former for small mi 
and m.2 irrespective of the values of a, g\ and g 2 ■ Explicit computations of 5N and SN C to second order in the field 
fluctuations are carried out in Appendix A. 



IV. CURVATURE PERTURBATION IN TWO-BRID INFLATION 



Let us evaluate the curvature perturbation explicitly. The quantities of our interest are the curvature perturbation 
spectrum Vs, the spectral index ns, tensor-to-scalar ratio r, and the non-Gaussianity /jv£ al - To evaluate these 
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quantities, we expand the SN formula (3.13) to second order in 8<f> for N(<f>) given in (3.12). Note that £7 must be 
expressed in terms of S(f> with second order accuracy, using Eq. (3.11). Details are given in Appendix A. The result is 



XAT 0</>i5icos7 + (5(/> 2 £ 2 sin7 g ( g £ {m 2 8<t>\ - mrffay 

oJy = : 1- — — 7 : + • • • . (4.1) 

migi cos 7 + m 2 g 2 sin 7 la \m\g\ cos 7 + m 2 g 2 sin 7) 13 

We assume that the scalar field fluctuations 6<fri and S(j) 2 are Gaussian, with the dispersion, 

Sab , (4.2) 



{5(j} A 5(pB) k = ( ^ 



H ,2 



where tk is the horizon-crossing time of the comoving wavenumber k, k = Ha. Then the curvature perturbation 
spectrum is given by 



-Ps{k) = ^Pn(k) 91 C ° S2 7 + 91 Sin2 7 



(2-7t) 3 (miji cos 7 + m 2 g 2 sin7) 2 \2ir J 

Using the fact that 3H 2 = V, the spectral index is easily calculated as 

d\nVn(k) 1 , 2 , 2n 



(4.3) 



= 1 + rflnfc v = 1 - (mf + mi) . (4.4) 



Thus our model predicts ns < 1. As for the tensor-to-scalar ratio, it is given by 

^ _ Vt_ _ g (miffi cos 7 + TO 2 g 2 sin 7) 2 
"Ps cos 2 7 + g\ sin 2 7 



(4.5) 



As we see from Eqs. (4.3), (4.4) and (4.5), our model seems to have a sufficient number of parameters that can be 
tuned to give the values of ns and r which are consistent with observations. This is in contrast to a single-field hybrid 
inflation model with V = Vq e m ^, for which the number of e- folds is given by 

iV singlc - , (4.6) 

which gives 

ns = 1 — m 2 , r = 8m 2 (single-field case) , (4.7) 

and, to first order in the slow-roll approximation, the non-Gaussianity is exactly zero. Comparing this single-field 
case with our two-field case, it is amazing that adding only one extra component to the inflaton field results in the 
huge variety of the output parameters. 

Now we evaluate the non-Gaussianity in our model. For convenience, we introduce the linear curvature perturbation 
TZl and the linear entropy perturbation S, 

<5</>i3icos7 + <S</> 2 g 2 sin7 Sfa g 2 sin 7 - 5<t> 2 gi cos 7 

l<-L = ; : , o = : . (4.8) 

m\g\ cos 7 + m 2 g 2 sin 7 m 2 gi cos 7 — mi g 2 sin 7 

For the Gaussian fluctuations 84>a given by Eq. (4.2), we see that they are orthogonal to each other: 

{H L -S)=0. (4.9) 

In terms of TZl and 5, the nonlinear SN in Eq. (4.1) is expressed as 

SN = K L + \f£T\K L + Sf , (4.10) 
5 



where 



.local = Sgiffl (m 2 gi cos7-TOig 2 sin7) 2 

NL 6cr(g 2 cos 2 7 + g 2 sin 2 7) 2 m\g\ cos 7 + m 2 g 2 sin 7 
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Since both mi<7iC0S7 and 771232 sin 7 are positive, we thus see that in this model /j^f 3,1 is always non-negative. In 
addition, we note that this non-Gaussianity comes essentially from the end of inflation, as seen from its explicit 
dependence on the parameters that define the end of inflation. This fact may be explicitly seen by separating 6N into 
two parts; one from the perturbation during inflation 5N* and the other from the end of inflation 5N e . This is done 
in Appendix A. 

We also see that the linear entropy perturbation S induces a curvature perturbation at second order. It may be 
worthwhile to note that if one finds a way to extract out observationally the second order term, it may be possible 
to detect this second order effect of the entropy perturbation. In this connection, we mention recent work by Enqvist 
ct al. [21] and by Chambers and Rajantie [10], in which they discussed a mechanism of generating non-Gaussianity 
at the end of inflation through parametric resonance when there exist two light scalar fields. Their mechanism 
corresponds effectively to the case where the second order term in (4.10) is dominated by S 2 . In this case, the 
curvature perturbation may be approximately expressed as 

6N = K L + f s S 2 . (4.12) 

Chambers and Rajantie [10] define a non-Gaussian parameter essentially by the ratio, 

(4-13) 

Then, under the assumption that TZl dominates over fsS 2 in SN, that is, /| (S 2 ) <C (^|), one finds 

fNL ~ f5 Wf- (4 ' 14) 

In our case, since f s ~ /]^ al and (S 2 ) ~ (U 2 L ), our .f™ wou id be small - (f^ 1 ) 3 (S 2 ) even for f° c L al ~ 100. 

To examine the viability of our model, let us compare our predictions with observations. Analysis of the WMAP 
5-year data [22, 23] gives the spectral index of the value, 

n s = 0.96 ±g;g!4 ; (4-15) 
and the upper limit on the tensor-to-scalar ratio, 

r < 0.43 (95% CL: WMAP alone) , 

r < 0.20 (95% CL: WMAP+BAO+SN) . (4.16) 

As for the non-Gaussian parameter /]^£ al , the situation is somewhat controversial. Yadav and Wandelt [24] finds 
non-vanishing non-Gaussianity in the 3-year WMAP data [25] as 

27 < /nT 1 < 147 (95% CL) , (4.17) 

while Komatsu et al. finds only a bound in the WMAP 5-year data [23] , 

-9 < j^2 al < 111 (95% CL) . (4.18) 

In any case, these results indicate that a positive value of /]^£ al is more favored over a negative value, if non-zero. 
This is perfectly consistent with our result, Eq. (4.11). 

To show that these observational data can be indeed reproduced in our model, let us consider a specific set of the 
model parameters. We set 

ml - 0.005, m\~ 0.035, (4.19) 

and assume that the inflationary trajectory satisfies 

mi cos 7 3> m 2 sin 7 , (4.20) 

that is, 7 <C 1. We also assume that the coupling constants g\ and gi are of the same order of magnitude. In this 
case, ns and r depend only on mi and m2, 

n s = 1 - (ml +m 2 2 ) - 0.96, r w 8m^ - 0.04 . (4.21) 
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On the other hand, the non-Gaussian parameter /]^£ al depends on the other model parameters. For simplicity, let us 
assume 

g\ = g\ = g 2 ■ (4.22) 

Then the curvature perturbation spectrum is approximately given by 

^4(i) 2 ~ 2 - 5xI °" (4 23> 

where the second equality is from the WMAP normalization [26] . Since we have fixed m\ already, this normalization 
determines H 2 , and hence a 4 as a function of A, 

# 2 ~40m^ s ~5xl(r 10 , a 2 = (12\H 2 ) 1/2 ~ A 1/2 x 1(T 4 . (4.24) 

This gives 

f local _ 5 ff m 2 _ An _9_ (a 9^ 

Thus, in particular, the non-Gaussian parameter /]^£ al can be large if A is very small. Note, however, that A cannot 
be too small. For a » H which is necessary for the field % to work as a water-fall field, we must have A 1 / 2 3> 1CP 6 
from Eqs. (4.24), hence A 1 / 4 > 1CT 3 . 



V. CONCLUSION 



We have investigated analytically the curvature perturbation from a model of multi-field hybrid inflation, which we 
call ll multi-brid inflation" . 

First we have considered a general condition for a model to be exactly soluble under the slow-roll approximation. 
Then we have presented a multi-brid inflation model in which the curvature perturbation is analytically computable 
to full nonlinear order. We have described a method to calculate the curvature perturbation analytically. Here we 
have noted that the coupling of the (multi-component) inflaton field to a water-fall field that terminates inflation 
plays a significant role in the determination of the curvature perturbation. 

Then as a specific example, we have focused on a two-brid inflation model, and derived an analytical expression for 
the curvature perturbation to second order in the field fluctuations expressed in terms of the model parameters. Using 
this expression, we have given analytical formulas for the spectrum of the curvature perturbation Vs, the spectral 
index ns, the tensor-to-scalar ratio r, and the non-Gaussian parameter /]^£ al . We have noted that adding only one 
extra component to the inflaton field increases substantially the degrees of freedom in the output parameters that arc 
observationally constrained. We have found that our two-brid model can explain the recent WMAP observations [22, 
23, 25]. 

One possibly subtle issue we have not discussed in this paper is the effect of isocurvature perturbations at or 
after reheating. When the universe is reheated by the water-fall field (and perhaps by the inflaton field as well), 
the abundance of created particles may depend on the initial values of the components of the inflaton field. This 
can be regarded as a two-field version of the modulated fluctuations [6, 7]. This would lead to additional power to 
the curvature perturbation spectrum and probably to additional non-Gaussianity. This is an interesting effect which 
deserves further study. But from the point of view of the model we discussed in this paper, this dependence must be 
weak enough for our model to be viable. This effect must be carefully examined when we actually attempt to carry 
out the model-construction. In fact, Barnaby and Cline found that the effect is significant and leads to non-trivial 
constraints on some models of hybrid inflation [27] . 

Recently a lot of efforts have been paid in the construction of stringy inflation models (see e.g. [28, 29] and references 
therein; for reviews on stringy cosmology, see e.g. [30]). Many of these models fall into a class called brane inflation [31] , 
in which a brane is attracted to another brane (or anti-brane) in higher dimensions and the distance between the 
branes play the role of the inflaton. In these models, inflation ends with a collision of the branes. Then the collided 
branes annihilate to heat up the universe. Thus, in the 4-dimensional effective picture, brane inflation is the same as 
hybrid inflation. Therefore it seems quite possible to construct a stringy model that gives rise to multi-brid inflation. 
This is currently under investigation. 
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APPENDIX A: CALCULATION OF SN TO SECOND ORDER 



Let us first evaluate the perturbation in 7 to second order. Setting Sj = £17 + 5 2 1, where S\j and S 2 j are of linear 
and second orders, respectively, we take the perturbation of Eq. (3.11) to second order, assuming Scfii and S(p 2 are of 
linear order. We have 



mi 



0<P2 

m 2 



sin 7 cos 7 



gimi g 2 m 2 



(617 + 627) - 



a I cos 7 



sin 7 



2 \5i™i 52«i 2 j 



if 



The linear part of the above equation determines £17. We find 

gig 2 m 2 S(pi — miS(j) 2 



5i7 : 

a g\m\ cos 7 + g 2 m 2 sin 7 

Then collecting the second order terms in Eq. (Al), we find 

1 gimx sin 7 - .92^2 cos 7 2 
02I = 77 ; : (Oil) 



(Al) 



(A2) 



2 gimi cos 7 + g 2 m 2 sin 7 

g 2 m 2 cos 7 



.9i.92 9i m i sin 



2cr 2 (gimi cos 7 + g 2 m 2 sin7) 3 



[m 2 84>i - rriiS4> 2 ) 2 . 



(A3) 



Now it is straightforward to calculate SN. Expanding eq. (3.12) to second order in the perturbation, and substituting 
Eqs. (A2) and (A3) into the resulting equation, we obtain 



SN : 



gi cos 7 + 8<j>2 ff2 sin 7 g\g\ (m 2 o>i - m 1 84> 2 )' 1 



toi<?i cos 7 + m 2 g 2 sin 7 2cr (toi^i cos 7 + m 2 g2 sin7) 3 



(A4) 



This is recapitulated in Eq. (4.1). 

It may be noted that as long as explicit computations of 8N are concerned, it is easier to go back to the original 
slow-roll equations (3.2). Integrating them are trivial. We immediately find 



01 



N . 



02 



& 2,/ 



N. 



mi mi 
where <j)\j and <f) 2 j are parameterized in terms of 7 as in Eq. (3.10), 



h.f = — cos 7, 
9i 



t> 2 f = — sin 7 . 

52 



(A5) 



(A6) 



One can then take the perturbation of either the first or the second equation in Eq. (A5) to second order, with 8<j)ij 
and 8(j) 2 .f given by 



HiJ = — 
9i L 



a 
91 



sin 7 [8\i + ^27) — cos 7 
cos 7 (<5i7 + S 2 j) — sin 7 



(til) 2 



2 

(<5i7) 2 



(A7) 



The result is the same whatever method one adopts. In fact this may be used to check the calculation. For example, 
taking the perturbation of the first equation in ( A5) and using the first line of Eq. (A7) , one finds 



8N 



8(f>\ 8<j>ij 



mi 



mi 



001 a . . a 
1 sin7oi7 H 



mi 771151 



migi 



sin 7 (527 + cos 7 



{hi? 



(A8) 
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Inserting Eqs. (A2) and (A3) into this equation, one can check that the result agrees with Eq. (A4). 

Now let us turn to the evaluation of the correction term after inflation, SN C . Taking the perturbation of N c in 
Eq. (3.15), we find 



cat 1 r r i f i l a I m\ . m 2 

oN c = - [midcpij + m20(p2j\ = —-t \ — sin 7 COS7 017 

4 4 \ g\ g 2 



mi . m 2 \ f mi m 2 . \ ($iy) 2 
sm 7 cos 7 I d 2 7 + cos 7 H sin 7 1 



(A9) 



cr 

4 L V 9i 92 J \9\ 92 

Inserting <5i7 and <5 2 7 given by Eqs. (A2) and (A3) into the above, we obtain 

XAT lTO lff2 sin7-m 25l cos7 ffitf! ( m i + m|)(m 2 <tyi - mi5<j> 2 ) 2 ,. . 

6N C =- ■ : — (m 2 o<pi - mx6<p2) ? ■ : r^— • (A10) 

4 mi 51 cos 7 + m 2 g 2 sin 7 8ct (migi cos 7 + m 2 g 2 sin7) ;1 

Comparing this with Eq. (A4), we see that the contribution of SN C to <JAT is suppressed by a factor quadratic in mi 
and/or m 2 . Hence, under the assumption that mi and m 2 are sufficiently smaller than unity, which is necessary for 
the slow-roll condition to hold, we may neglect the correction term SN C . 

Finally, let us evaluate SN up to a surface of constant potential V = const. = V* during inflation, which we denote 
by SN*, and identify the contribution from the end of inflation, which we denote by SN e . On super horizon scales, 
and uder the slow-roll approximation, SN* is equal to the curvature perturbation on the uniform density slice, which 
is usually denoted by (. 

The slice V = V* is given by 

mi^i, +m 2 2 ,» = ln(K/Vb) = C. (All) 
One can parametrize 4>a,* {A — 1, 2) with a parameter t as 

1 ,„ 1 



(C + 2t), <i>2,* = 7 , — (C-2t). (A12) 



2mi ' 2m 2 

This gives 

(501 * = — Jt , 5<j)2 * = — — <ft . (A13) 
mi ' m 2 

Thus taking the perturbation of the solutions (A5), with <j)A,f replaced by <f>A,* and N hy N*, we obtain two expressions 
for SN*, 



H 1 _ 


Hi,* 


_6fa 


St 


mi 


mi 


mi 


m\ 


5<t> 2 


5<p2,* 


_ 5<j)2 


St 

+ — 

mi 


m 2 


m 2 


m 2 



From these equations, we find 



SN* = = (A14) 



St = ^ im2 2 {m 2 5(f)i — miS<j} 2 ) . (A15) 



m\ + m 2 

Inserting this to either one of Eqs. (A14), we obtain 

SN* = mi ^ + ^ ■ (A16) 
m\ + mi 

As clear from this expression, in our model, there is no non-Gaussianity from the evolution during inflation. This is 
in agreement with the discussion based on the equations in the qa space given near the end of Section III. 

Subtracting SN* from SN given in Eq. (A4), we obtain the contribution from the end of inflation SN e . The result 

is 

5N e = SN-SN* 

(m 2 g 1 cos'j ~ mig 2 sm-f)(m 2 5(j) 1 - mi5(f>2) g\g\ (m 2 S(f>i - miS(f> 2 ) 2 



+ v z yi 1 yz : .. . (A17) 



(m\ + m 2 )(mi<7i cos 7 + m 2 <? 2 sin 7) 2a {m\gi cos 7 + m 2 g 2 sin 7) 
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